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1.  Introduction 


,'~^n  this  paper  we~est*bi  Lsh  some  refined  estimates  for  the  approximation 
of  the  eigenvalues  and  eigenvectors  of  selfadjoVnt  eigenvalue  problems  by 


/ 


finite  element  or,  more  generally,  Galerkin  methods.  Suppose  \  is  an 


eigenvalue  of  multiplicity  q  of  a  selfadjoint  problem  and  let  M(A)  denote 


the  space  of  eigenvectors  corresponding  to  A. ~\  Denote  by  fl*||g  the  energy 


norm  for  the  problem.  Let  ^^^0<h  ^  of  finite  dimensional 


approximation  spaces  employed  in  the  Galerkin  method.  A  will  be  approximated 
from  above  by  q  of  the  Galerkin  approximate  eigenvalues: 


ASA.  .  S***S  A. 

h,  1  h,  q 


■  J  f 

6 


>>,/  "  So 


<'  F 


A  a  A  ,  •  •  • ,  A  . 
h,  1  h,q 


Let  u,  with  ||ujg  =  1,  denote  an  eigenvector  corresponding  to  A,  and  let 

Uj^  q*  with  k||B  *  1,  denote  the  Galerkin  eigenvectors 

corresponding  to  A  ,  •••, A.  ,  respectively. 

n,  i  n, q 


It  is  well-known  that 


(1.1) 


A,  ,-A 
h,  k 


S  C  sup  inf  |u-*||g  ,  k  *  l,***,q, 


ueM(A)  %eS 


Mb"1 


and  that  there  is  a  u^  =  u^Ch)  €  M(A),  with  flu^g  =  1.  such  that 


(1.2) 


liuh  k"  "Jr  5  C  sup  inf  ||u-z||  ,  k 

ueM(A)  *eS. 

Mb'1 


=  1, ***,q. 


In  [7,8]  Chatelin  proved  the  following  refinements  of  (1.1)  and  (1.2): 
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(1.3a)  |u  -  E^ujg  =  i{a)  Inf  |tt-z|B  V  u  €  MU) 


xeS. 


(b) 


(1.3b)  IVk  '  Euh.k*B  =  rh  ‘f  k-l.  —  .q 


and 

(1.3c) 


xeSv 


l<\t  -  X)/X«B  *  ^C)  >f  iE"h.k-*lr  k  = 


where  E  denotes  the  orthogonal  projection  of  the  energy  space  onto  M(X) 
and  E^  the  orthogonal  projection  onto  span{u^  !•***>%  and  where 

f  i) 

r,  — »  1  as  h  — >  0,  for  l  -  a,b,c. 

h 

The  purpose  of  this  paper  is  twofold.  The  first  purpose  is  to  establish 
an  estimate  for  |r^-l|.  We  show  that 


(1.4) 


|r*£)-l|  £  dr)2(h). 


where  i?(h)  is  a  certain  measure  of  the  approximability  property  of  {S^}; 
for  the  definition  of  v  see  Section  3.  This  is  done  in  Section  4. 

In  (3)  the  authors  established  the  estimate 

(1.5)  \  -  X  =  C  inf  inf  ||u-x||2  , 

n>1  ueM(X)  *€S, 

UIIB=1 

which  is  an  improvement  over  (1.1)  and  (1.3c)  in  the  case  of  a  multiple 
eigenvalue.  [3]  also  contains  estimates  for  X^  ^-X, k=2, • • • , q,  and  for 
||u^  k  *  1,  which  are  improvements  of  (1.1)  and  (1.3c)  and  of 

(1.2)  and  (1.3a,b),  respectively.  The  second  purpose  of  the  paper  is  to 
present  a  simplified  proof  and  an  extension  of  the  results  in  [3],  This  is 
done  in  Section  5. 
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In  Section  2  we  give  a  precise  statement  of  class  of  eigenvalue  problems 
and  approximation  methods  we  will  consider.  Section  3  contains  some  back¬ 
ground  information. 

The  2nC*  author  would  like  to  thank  Professor  Hans  F.  Weinberger  for 
several  helpful  discussions  on  the  topics  in  this  paper. 


Setting  for  the  Problem 


Suppose  H  is  a  real  Hilbert  space  with  inner  product  (•,•)  and  norm 


U *  || ,  respectively,  and  suppose  we  are  given  two  symmetric  bilinear  forms 
B(u,v)  and  D(u,v)  on  H  x  H.  B(u, v)  is  assumed  to  satisfy 


(2.1)  |B(u, v)  |  £  Cjuf  flvfl,  V  u.v  €  H 
and 

(2.2)  C0||ufl2  *  B(u,u),  Vu  e  H,  with  C2  >  0. 

It  follows  from  (2.1)  and  (2.2)  that  (u,v)g  *  B(u,  v)  and  ||u||B  *  B(u,u)1/2 
are  equivalent  to  (u.v)  and  ||u||,  respectively.  Regarding  D  we  assume 


(2.3) 
and  that 

(2.4) 


is  compact  with  respect  to  ||*||,  i.e.,  from  any  sequence  which  is  bounded  in 
||*  j| ,  one  can  extract  a  subsequence  which  is  Cauchy  in  ||*||q-  For  the 
remainder  of  this  paper  we  will  use  B(u,v)  and  ||*||g  as  the  inner  product 
and  norm  on  H  and  denote  this  space  by  Hg. 

We  then  consider  the  variational ly  formulated,  selfadjoint  eigenvalue 
problem. 


3 


(2.5) 


Seek  A  (real)  and  0  *  u  €  Hg  satisfying 
B(u,v)  =  AD(u,v),  V  v  €  Hg  . 


Under  the  assumptions  we  have  made,  (2.5)  has  a  sequence  of  eigenvalues 

0<A1  S  Ag  S  •••  /"  +  oo 

and  corresponding  eigenvectors 

U1 ’ u2’ *  *  * ’ 

which  can  be  chosen  to  satisfy 


(2.6) 


Btu^Uj)  =  A^Dtu^.Uj)  =  5ij,i,j=l,2, 


The  eigenvalues  auid  eigenvectors  satisfy  the  following  well-known  variational 
principles: 


(2.7) 


and 


min 


B(u,u)  .  B(W 
D(u,u)  DTuj^Tu^y 


B(u, u.  )=0 

i=l ,2, • • • , k-1 


,  k  = 


1.2, 


(the  minimum  principle) 


(2.8) 


,  B(u,u) 

*  min  max  h7 — — r  =  max 

\cHb  ueVk  U'U  ueUk*sp(ur 

dim  V.  =k 
k 


’V 


B(u,  u) 
D(u, u) ’ 


k=l , 2, •  •  • 


(the  minimum-maximum  principle). 


For  any  A^  we  let 


(2.9)  M  =  M(A^)  *  {u:  u  is  an  eigenvector  of  (2.5)  corresponding  to  Ak> . 


We  shall  be  Interested  in  approximating  the  eigenpairs  of  (2.5)  by 


finite  element  or,  more  generally,  Galerkin  methods.  Toward  this  end  we 


suppose  we  are  given  a  (one  parameter)  family 


hi 


of  finite  dimensional 


0<h£l 


subspaces  c  Hg  and  we  consider  the  eigenvalue  problem, 


(2. 10) 


Seek  Aureal),  0  *  u^  e  satisfying 

BCu^v)  =  Xj^Dtv^.v),  V  v  €  Sh  . 


The  eigenpairs  (A^.u^)  °**  (2.10)  are  then  viewed  as  approximations  to  the 

eigenpairs  (A.u)  of  (2.5).  (2.10)  is  called  the  Galerkin  method  determined 

by  the  subspaces  {S^>  for  the  approximation  of  the  eigenvalues  and  eigen¬ 
vectors  of  (2.5).  We  will  also  sometimes  refer  to  problem  (2.10)  as  the 
Galerkin  approximation  of  problem  (2.5).  (2.10)  has  a  sequence  of  eigenvalues 

0<Vl“h.2  S"'SVk’  N-dl"  Sh’ 


and  corresponding  eigenvectors 


VrVa’^-VN 


which  can  be  chosen  to  satisfy 


(2. 11) 


B(u 


u.rVj 


>  =  ViD(Vi'Vj)  ■ 


5^,  i,  j=l,  •• 


N 


The  (A^  j,uh  are  referred  to  as  the  approximate  eigenpairs,  while 
(X.,u. )  are  referred  to  as  the  exact  eigenpairs  of  (2.5).  Maximum  and 
minimum-maximum  principles  analogous  to  (2.7)  and  (2.8)  hold  for  problem 
(2.10);  they  are  obtained  from  (2.7)  and  (2.8)  by  replacing  by  Sh  and 
letting  k=l,***,N.  We  will  refer  to  them  by  (2.7^)  and  (2.8h),  respective¬ 
ly.  Using  (2.7)  and  (2.8)  together  with  (2.7*1)  and  (2.8^)  we  see  immediately 
that 

(2.12)  Xk  "  Xh  k’  k  =  =  dim  Sh' 

We  will  assume  that  the  family  {S^>  satisfies  the  approximabi 1 ity 
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assumption 

(2.13)  c  (h)  =  ||u§  1  inf  ||u-;t||g  — *  0  as  h  — »  0,  for  each  u  €  Hg. 

B  ,«Sh 

It  follows  from  (2.7),  (2.8).  (2.7h),  (2.8h),  and  (2.13)  that 

(2.  14)  A,,  — »  A,  as  h  — *  0,  for  each  k. 

h,  k  k 


Finally  we  introduce 


“j  '  ^  uj 


the  exact  eigenvectors  normalized  in  ||*||q>  and 


\,j  ^h,  j  “h,  j  ’ 

the  approximate  eigenvectors  normalized  in  H'flo- 

Throughout  the  paper,  the  specific  eigenfunctions  satisfying  (2.6) 
((2.11))  will  be  denoted  by  u^u^  j).  Thus  the  u^(u^  j)  are  normalized  in 
|| *  || g ;  u^fi^  denotes  the  same  eigenvectors,  renormalized  in  ||'||D-  When  we 
denote  an  eigenpair  by  (A,u)  we  will  not  assume  any  particular  normalization 
on  u.  C,Cj,d,  and  d^  will  denote  generic  constants. 
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3. 


In  this  section  we  present  several  preliminary  results  that  will  be  used 
in  the  sequel.  For  further  information  on  eigenvalue  problems  we  refer  the 
reader  to  [4,83. 

a)  An  Identity  Relating  the  Eigenvalue  and  Eigenvector  Errors 

Here  we  present  an  identity  that  relates  the  errors  in  eigenvalue  and 
eigenvector  approximation. 

Lemma  3.1.  Suppose  (A,u)  is  an  eigenpair  of  (2.5),  suppose  w  is  any  vector 
in  Hg  with  8wj|D  -  1.  and  let  A'  =  B(w, w).  Then 

(3.1)  A' -A  =  J|w-u||g  -  A|w-u|*  . 

Proof.  By  an  easy  calculation, 

(3.2)  I*-ul§  -  =  ||u||  -  2B(  w,  u)  <■  |u|^ 

-  M“Id  *  2*  °<“.“>  ♦  *Md  ■ 

Now 

B(v, u)  =  AD( v, u)  V  v  e  Hg, 

from  which  we  get 

(3.3)  B(w,u)  =  AD(w,u) 
and 

(3.4)  |u||g  =  B(u,  u)  =  AD(u,  u)  =  A|ju||p  . 

The  result  follows  from  (3.2)-(3.4)  and  the  relations  A'  =  ||w||g  and 

'  *  Md  • 

b)  The  Operators  T  and  Tk 
Let 
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Hp  =  the  completion  of  Hg  with  respect 


Hg  is  a  Hilbert  space  with  inner  product  D  and,  since  i*!g  is  assumed  to 
be  compact  with  respect  to  fl'flg.  Hg  is  compactly  imbedded  in  Hg. 
(Alternatively,  we  could  have  assumed  Hg  c  Hg,  compactly,  and  let  D(u,v)  be 
the  inner  product  on  Hg. ) 

From  Hg  and  Hg  construct  the  “negative  space"  Hg  =  H^  ,  with  norm 
|| *  |Lg-  Then  Hg  c  H_g  compactly,  and  for  v  €  Hg,  D(u,v)  has  a  continuous 
extension  to  u  e  H  g  so  that  D(u, v)  is  continuous  on  Hg  x  Hg.  For 


u  e  H_  ,  ||u||_B  =  sup  1 — 1  .  For  a  complete  discussion  of  this  construc- 

■"  V€Hb  V  1 


II B 


tion  we  refer  to  (5,  pp. 31-39]. 

Next  we  introduce  the  operators  T,T^:  H  g  — »  Hg  defined  by 


(3.5) 


Tf  €  Hg 

[  B(Tf , v)  =  D(f , v) ,  V  v  e  Hg, 


(3.6) 


<  V  e  sh 

B(Thf.v)  •  D(f.v),  »  V  €  Sh- 


T  and  are  the  solution  and  approximate  solution  operators  for  the 

"boundary  value"  problem  corresponding  to  the  eigenvalue  problem  (2.5).  It 
follows  Immediately  from  (2.1),  (2.2),  and  the  fact  that  D(f,v)  is  continuous 
on  H_gX  Hg  that  T  and  are  bounded  from  Hg  to  Hg.  Since  Hg  is 

compactly  imbedded  in  Hg  and  Hg  is  compactly  imbedded  in  H_g,  T  is 
compact  from  Hg  to  Hg,  from  Hg  to  Hg,  and  from  H_g  to  H_g.  is,  of 

course,  also  compact  on  Hg,  Hg,  and  H_g.  It  is  easily  seen  that  T 
and  are  selfadjoint  on  Hg  and  that  T  is  selfadjoint  and  positive 

definite  on  Hg  (with  respect  to  B(u,v)).  It  is  immediate  that  T  has 

eigenvalues 
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and  eigenvectors 


»*1  =  \l  *  ^2  =  X~2  ***  \  0 


U  11  •  •  • 

1’  2’ 


and  that  has  eigenvalues 


=  X:‘-  *  •••  *  «h,N-  •  N=s  dl”  Sh’ 


^h, 1  Xh, 1 


and  eigenvectors 


“h,  1’  *  *  *  ’  Uh, N‘ 

Let  Ph  be  the  orthogonal  projection  of  Hg  onto  S^;  then  from  (3.6)  we 
see  that 


Let 


(3.7) 


and 


rj(h) 


II  ( I — Ph )  I II Hg-^Hg  *  HT 


sup  inf  ||Tg-z|L 
«6HD  *eSh 

b«iid=i 


(3.8) 


u(h)  =  |(I-Ph)T|| 


P-Th| 


V«B 


sup  inf  ||Tg-*||  . 

*€HB  *eSh 


Several  of  the  results  in  Sections  4  and  5  are  stated  in  terms  of  the 
qualities  of  y  and  v.  We  now  present  some  properties  of  y  and  v. 


Lemma  3.2.  There  are  positive  constants  and  such  that 

(3.9)  Cji>(h)  £  tj(h)  <  C2  vV(h)  . 

Proof.  Since  ||u||g  £  C |J u J| ^  V  u  e  Hg  we  have  p(h)  S  Cr)(h),  which  is  the 
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- 1 

first  inequality  in  (3.9)  with  =  C  .  Now  consider  the  second  inequality 
in  (3.9).  From  (3.5)  and  (3.6)  we  have 

imb  *  M-b  •  nyis s  iifLB 

and  hence 


(3. 10) 

and  from  (3.8)  we  have 


<  2 


(3. 11) 


i>(h) . 


We  now  note  that  H  and  H0  are  connected  by  a  scale  of  Hilbert  spaces. 

~D  D 

It  thus  follows  from  (3.10),  (3.11),  and  a  result  on  interpolation  of  linear 
operators  [5,  pp.  240-242]  that 

"<h)  *  »T-Th»HD-^B  2  C2'/2  =  Cl"h)1/2’ 

which  is  the  second  inequality  in  (3.9).  □ 


Lemma  3.3. 

(3.12)  lim  Tj(h)  =  lim  r(h)  =  0  . 

h-*0  h-*0 


Proof :  Because  of  Lemma  3.2  it  is  sufficient  to  show  that  lim  t>(h)  =  0. 

h-£> 

(2.13)  implies  that  P  — >  I  pointwise  on  HD  (in  fact,  (2.13)  is  equivalent 

n  d 

to  this  result).  Since  T:Hg  — »  Hg  is  compact,  T(geHg:  ||g||g  =  1}  is 

relatively  compact  in  H^,  and  lim  i>(h)  =  0  follows  from  the  standard  result 

B  h-?0 

that  a  family  of  operators  that  converges  pointwise  on  a  space  converges 
uniformly  on  a  relatively  compact  subset.  □ 
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2  2 

From  Lemma  3.2  we  have  .  tj  £  OCv).  It  may  happen  that  tj  =  o(i>).  This 
is  shown  by  the  following  example. 

Example 

Let 

Hg  =  Hq(0,  1 ) , 


and 


B(u, v) 


u' v'dx, 


D(u,v)  =  J  u  v  dx, 

0 

where  0  <  a  <  a(x)  <  £  <  «.  ( (0,1)  is  the  l  order  Sobolev  space  and 
Hq ( 0 , 1 )  =  {u  e  H1 (0, 1 ) :  u(0)  =  u(l)  =  0> ) .  For  f  e  L2(0,1),  u  =  Tf  is  the 
solution  of 


’  -(a(x)u' ) '  =  f  ( x) ,  0  <  x  <  1 
u(0)  =  u(l)  =  0. 


First  suppose  =  the  space  of  C^,  piecewise  linear  functions  with 

mesh  size  h  that  vanish  at  0  and  1  and  suppose  a(x)  is  smooth.  Then 

2 

we  easily  see  that  ij(h)  ~  h  and  v(h)  ~  h,  so  that  tj  =  o(v).  Next 
suppose  =  the  space  of  C*"*,  piecewise  quadratic  functions  vanishing 

2 

at  0  and  1.  If  a(x)  is  smooth  we  see  that  rj(h)  ~  h  and  r(h)  ~  h  ,  so 
2 

tj  ~  v.  If,  on  the  other  hand,  a(x)  is  rough,  specifically  if  a(x)  is 

such  that  g  €  Hp  =  LA0,1)  implies  u  =  Tf  €  H^(0,1),  but  g  e  Hg  =  Hq (0,1) 

does  not,  in  general,  imply  u  6  Ha(0,l)  for  a  >  2,  then  tj  ~  h  and  v  ~  h, 

2  (  ^ 

SO  TJ  =  O  (  V )  . 

From  (2.13)  we  have 
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|U-Ph)u|B  -  eu(h>  |u|B  -.0,  V  u  6  Hg. 

The  usual  dually  argument  (cf.  Aubin  (1],  Nitsche  [10],  and  Oganesjan- 
Rukhovets  [11])  shows  that  ^ ( I— S  Ctj( h)  |  ( I-P^Juflg  and 
|l(I-Ph)u|  g  £  Cu(h)  ||(  I-P^)u|jg.  For  the  sake  of  completeness  we  Include  proofs 
of  these  results. 


l(I-Ph)u|D  s  „(h)|(I-Ph)u|B.  v  u  6  Hg 

l(I*Ph)ul-B  £  »  u  <E  Hg. 


Lemma  3.4. 

(3. 13a) 
and 

(3. 13b) 


=  sup  |B((I-P  )u,Tg)| 
«eHD 

»s||D=1 
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which  is  (3.13a).  Similarly, 


S  sup  inf  j|Tg-%||B  |(I-Ph)u| 
g€HD  *eSh 

<  r,(h)|(I-Ph)uIB. 


fl(I-ph)uLB  =  SUP  |D((I-Ph)u.g)| 
geHB 

=  sup  inf  |Tg-zIB  ||(I-Ph)u|| 
g€HB  *eSh 

WB-1 

"  y(h)||(I-Ph)uj|B. 

which  is  (3.13b).  a 


c)  Preliminary  Eigenvector  Estimates 

For  i  a  1,2, •••  let  be  the  lowest  index  of  the  it^1  distinct 

eigenvalue  of  (2.5)  and  suppose  X,  has  multiplicity  q. .  Let  E  =  E(X,  ) 

kj  i  kj 

be  the  orthogonal  projection  of  Hg  onto  M(X^  )  and  let  E^  =  E^X^  )  be 
the  orthogonal  projection  of  H_  onto 


(3. 15)  =  M^(X^  )  =  the  space  of  eigenvectors  of  (2. 10) 

corresponding  to  Xh, ., j=0, • • • ,qi~l. 


Lemma  3 . 5  There  is  a  constant  such  that 

(3.16a)  lu-Eh<*klHB  £  Cll(I-Ph)ulB  ■  VueM,\»- 


(3. 16b)  l'J'Eh,Xk1)ulD  S  .  VueM(Xki), 


13 


and 


(3.  16c)  lu-Eh< Vul-B  *  Cl*(I'Ph)Ul-B  ’  V  U  *  MUk,!’ 

Proof.  Suppose  the  spaces  Hg.Hg,  and  H_g,  the  bilinear  forms  B  and  D, 

and  the  operators  T,T^, E,  and  E^  have  been  complexified  in  the  usual 

manner.  Let  be  a  circle  in  the  complex  plane  centered  at  ^  -  A^1 

enclosing  no  other  eigenvalues  of  T.  Then  for  h  sufficiently  small, 

H,  .  =  A,1,  ,  •••,4.  ,  ,  =  X.1.  ,,  but  no  other  eigenvalues  of  T,  , 

Hi,  kt  h.kj  Hi.kt+q,-!  h,ki+qi~l  °  h 


are  contained  in  T,  ,  and 

M 


(3. 17a) 

and 

(3. 17b) 


E(A. 


=  2tt  i  L  ( 

lct 


(z-T)"1  dz 


V\)msrL  (z-vldz 

*1 


These  are  the  usual  formulas  for  the  spectral  projections  associated  with  T 
and  and  Th  and  ^  ^  respectively  (cf.  [9,  Section  XI. 9]). 

Consider  now  the  proof  of  (3.16a).  Using  (3.17)  we  have 


(3. 18) 


»u-Eh(V,ulB  “  ■tE(V)  -Eh,V",lB 
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S  I2n  rad  (1^)1  sup  ^  |(T-Th)u| 


zerv 

0<h 


h  “B 

13  ^-FidTf~y 

*1 


v-1. 


•  pk)  _sup  l^-Th)"‘lHB^B  l(I-Ph,ulB’  V  u  6  MUk,>- 


zer, 

0<h 


*i 


In  the  last  inequality  we  used  the  relation  (T-T.  )u  =  (I-P.  )Tu 

h  h 


%(I'Ph)u- 


Now  IT-TJ 


0  implies 


C1  ’  »V  l(Z'Th)'1|HB-»HB  < 


0<h 


so  we  have  established  (3.16a). 

Now  consider  the  proof  of  (3. 16b).  The  above  analysis  is  relative  to 
the  space  Hg  (the  integrals  in  (3.17)  converge  in  the  operator  norm  on  Hg 
and  Since  T  and  Th  can  also  be  considered  on  Hg 

|| T— Th || ^  can  aPp1^  the  same  argument  in  Hg.  Note  that  the 

formulas  (3.17)  will  now  define  projections  on  Hg  which  are  extensions  to 
Hg  of  E  and  E^.  We  thus  obtain  (cf.  (3.18)) 

ki 


and 


lU'EhWk,)UlD  5  "it,  z=“P  l<Z-Th,_l(HD^HD  l(I-Ph,U#D-  V  “  €M(V’ 


h>0 

which  is  (3.16b). 

The  proof  of  (3.16c)  is  similar. 


Remark  3.1.  It  is  essential  in  the  Lemma  3.5  that  h  is  sufficiently  small, 
meaning  small  in  comparison  with  the  gap  between  A  and  A.  ,  A 

Ki  kj-1  kj+1 

If  this  gap  is  small  then  it  caui  happen  that  the  approximation  eigenfunction 
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“h,  ki 


associated  with  X^  ^  could  be  close  to  or  u^^. 

Lemma  3.5  is  an  eigenvector  estimate  since  it  provides  an  estimate  for 

u(an  exact  eigenvector)  -  E^u(a  linear  combination  of 

approximate  eigenvectors). 


We  note  that  (2.13)  and  (3.16)  imply  that  )  is 

one-to-one  and  onto  for  h  sufficiently  small. 

We  next  prove  a  refinement  of  (3.16a)  due  to  Chatelin  [7,8].  (3.16a) 

shows  that 

E.  (X  )uj 

-*0(i>>  Vu€M(V- 

Chatelin  showed  that 


lu-Ehuk,  )u1b  , 


as  h  — »  0(see  (1.3a)); 


her  argument,  in  fact,  establishes 


Lemma  3.6 

(3. 19) 


(Chatelin) 

1  S 


There  is  a  constant  d. 


l 


lU~EhUk,  )U«B 


1  +  d^fh). 


such  that 


V  u  €  M(X  ), 

M 


where  u(h)  is  defined  in  (3.8). 


Proof.  For  the  sake  of  completeness  and  to  establish  the  form  of  the  bound  in 

the  second  inequality  in  (3.19)  we  present  a  proof  of  this  result. 

Let  T  =  P  T  P.  =  T.  P  .  Note  that  T,  and  T,  have  the  same 
h  h  h  h  h  h  h 

nonzero  eigenvalues,  that  E^tX^  )  commutes  with  T^,  and  that  Th  is 
selfadjolnt  with  respect  to  B.  For  u  €  M(X.  ), 
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(Th  -  "h.k.'V  -  phT!Ph-I)u  *  "W k.’V 


and  hence,  since  )  commutes  with  , 

i3-2o>  <fh  -  w’-vv’v*  *  ii-Eh(\1)iphTiph-i,u 

*  (pkl-^lcl,n-Eh,\1,,Phu- 


Let  Q  be  the  orthogonal  projection  of  Hg  onto  the  null  space 

of  Th-  Then,  any  z  e  SlI-E^fA^  )),  the  range  of  I-E^U^),  can  be  written 


as 


z  = 


N 

Z  B(z’uh.i)uh,j  *  Qz- 

t- 1 

l* ki>  ••*,lci+qi-l 


Here  we  have  used  the  orthogonal  decomposition 


Thus 


Hg  =  3J(Th)  ©  AXTh) 


=  x(t.  i  ©  m. ) 

n  n 


3  span  \,1'"-'VnI  •  *(V 

»  span  <Vr'"'VN>  ®  *<fh>- 


N 

(V\ k,)2  *  Z  Bl!’\,!l(Vrt,k,)\!  -"h.k,132’ 
1=1 

•••,ki+qi-l 


and  hence 


(3-21)  I'V'W^b  *  Z  lB(2-V«1l2K.t-Vk,l: 


1=1 

t*  kiP ••♦,ki+qi-l 
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*  l^.k.l2  Mb 

2  Bin  {|Vj  -  ^  J2,  J=1,--,N.J  .  V"’ W1’  K.k,!2} 

r  N  i 

X-  £  (BCz.i^ ^)|2  +  jQz|||  . 

1=1 

t*  kj, • * • , kj+q^-1 

» 

”ln  {K.k,-i  -  Vk,l2'K.k1.i  -  ^,k,l2  •K.kJ2}  »z*l-  1  2  2 
•  min  {K.k2  -  Vk.l2'  IVk.l2}  Mb'  1  =  »• 

Since  j  — »  (cf.  (2.14))  for  each  J  as  h  — »  0, 

■ 

"in  {lMh,k,-l"Mh,kJ  '  ^,k,*l  *  ’“h.k,  I  '  l“h,k1l  }'  122 

.  "ln  {K.ka-  *v.kj2 '  K.kil2}  '  1  =  1 

mln  {i‘^,-r'vll2'  K,.i  -  "k.l2  K,l2}  •  1 2  2 
"in  {Kr  “k,i2  •  K,l2}  •  1  * 1 

» 

2 

=  5^ ,  as  h  — >  0, 

from  (3.21)  we  get 

(3.22)  I'VVk.'^B  2  S1  Mb  '  V  b  a  *U-EhUki))  and 

V  small  h, 
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where  5^  >  0  and  depends  only  on  the  gap  between  and  fi^  +  1 

Combining  (3.20),  (3.22),  and  the  fact  that  I-E^(A^  )  and  are 

orthogonal  projections  we  have 


•  “-W  Vie  5  Cl'^'V  phT(ph-I)u 


♦  l%A.k,)(I‘Eh,V1  Vb 

s  *~i'  {»tiVi)2uIb  * 

l*Wk,l  l(I-Eh(V>  VIb}- 


from  which  we  get 


(3.23) 


l(I-EhtV»  Ph“«B  S  dllT(Ph-I,U«HB^B  I(Ph-IiU«B 

'  dtIIPh-I)TlHB^Bl(Ph-I)ullB- 


In  the  last  equality  we  used  the  fact  that  (P  -I)  and  T  are  selfadjoint 

h 

and  that  the  norm  of  an  operator  and  its  adjoint  are  equal. 

(3.23)  implies 

|l(I-Eh(Vph)u#B  ■  «(i-VuIIb|  5  »(I-Eh(V)phulB 


S  dil(Ph-I)TlHB^HB  »<Ph-I)UlB 


and  hence 


(3.24) 


J(I-Eh(\,)Ph,U#B  , 

i<ph-nuiiB 


5  dii(ph-I,Ti 


W 


We  easily  see  that 
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<(1-Ph,U»B  S  ■(I-Eh(\,,)U«B  S  l(I-EhUk,)Ph)UllB 


and  thus 


(3.25) 


.  ,  l(I-Eh(V)U«B.  «lI-EhUk,)V“«B 
'  I(Ph-I,ulB  "  «‘I-VuIb 


Combining  (3.24)  and  (3.25)  we  have 


SU-tyx  ))u|b 

0  5  -irph-i>u||B  - 1 


«(I-Eh(Xkl,Ph>u«B  , 

u  /  t  \  n  ~  ~  1 


«(Ph-I)ullB 


S  dll(Ph-IITlHB^HB  ■  V“«M(V- 


Recalling  that  |j(P  -I)T||h  „  =  u(h)  we  obtain  the  desired  result.  a 

n  “bhb 

Remarks  3.2.  (3.19)  should  be  compared  with  (4.21),  which  provides  a  stronger 

estimate  for  certain  special  u’ s  in  M(X,  ). 

M 

Lemmas  3.5  and  3.6  show  that  starting  from  any  exact  eigenvector  u  we 
can  construct  )u>  a  linear  combination  of  approximate  eigenvectors 

that  is  close  to  u.  One  can  also  start  with  an  approximate  eigenvector  and 
construct  a  close  exact  eigenvector.  We  present  another  result  of  Chatelin 
[7.8];  (see  (1.3b)). 


Lemma  3.7  (Chatelin).  There  is  a  constant  dj  such 


(3.26) 


1  5 


K.j  ~  , 

V^VVj-^VYjC 


j  =  kj.--.kj  +  q.  -  1. 
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Proof.  Observing  that  E(A^  )P^  =  (E^CX^^ J-E^X^ ))P^  +  E(X^^ ) C I-Ph) . 

we  obtain 

lE(\)-Eh(Ak1)PhlHB^B  £  *  «E<Xk,)!I-Ph)«HB^B' 

We  easily  see  that 


»EUk1)lI-Ph,»HB^B  *  l(I-Ph)E(VllVHB 

'  S“P  lU-Ph,EUk,)ulB 

ueHB 

Hb"1 


-  ^  i(i-ph,TEuk,)uiiB 

MB=‘ 

£  X.  u(h) 

M 


and  by  a  slight  modification  of  estimate  (3.18)  we  have 


Thus 

(3.27) 

Next  note  that 


lEUk,>-VXk1,«HB^B  S  Cp(h)- 

lE(\1,-EhUk,,PhlHB^HB  £  C‘'(h)- 


{'-'Wh  '  E(V’}  'Vj-e[V\j> 


*  (Eh(Vph-I1E(VVj' 


Hence,  using  (3.24)  and  (3.27),  we  have 

K,J  -  E<VVj'b  S  »{'-IEhUk1)Ph-EUk1)1}  1«Hb^b 
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X 


iIEh'Vph-I1E(Vv.j»B 

,  l(EhUk,)Ph-I1E(\,>VjlB 

„  "*<ll»ll(Ph-I)E(V%,jlB 

1-C» 

which  implies  the  second  inequality  in  (3.26).  The  first  is  immediate.  □ 


d)  Relation  between  Eigenvector  Error  in  H •  ||g_, _ Jj •  ||^,  and  || •  J|  . 

In  subsection  3.b)  we  noted  that  j(I-P^)u|{p  <  t>C h )  ||(I-P^)u||g  and 

||( I-P^)u||_B  £  u(h) |J( I-P^Jujg.  In  this  subsection  we  establish  similar 
results  for  the  eigenvector  error. 

For  i  =  1,2, •••  and  J  *  k.,*‘*,k,  +  q,-l,  let  u^  e  M(A  )  satisfy 

111  J  Ki 

E^tA^)  u|j  =  u^  j  .  We  know  from  the  discussion  in  Subsection  3.c)  that 
Uj  exists  and  is  unique  for  h  small.  From  (3.13a)  and  (3.16b)  we  have 

■  vA  ■  -  W^Id 


or 

(3.28a) 


s  C,  u(h)|(I-Ph)Uj|B 
S  Ci  ,,(h)IV"h.jlB  ’ 


Huj  '  "h.J*B 


S  CiD(h). 


It  follows  immediately  (by  scaling)  that 
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(3.28b) 


fluj  -  \  ,8d 

— £ - *  C.n(h)  , 

1“  -  Yjb 


where  u^  e  MU^  )  satisfies  E^U^JUj  =  j.  (Recall  that  ||u^  h||D  =  1 
and  ||Uj  h|JB  =  1.)  Similarly,  from  (3.13b)  and  (3.16c)  we  get 


(3.29a) 


Eli  -  Vll-B 

1=5  -  vjIb  1 


S  C, y(h) 


and 


(3.29b) 


fluj  *  %  jl-B 

— £ - 5  C,  w(h) . 


By  Lemma  3. 1  we  know  that 


(3'30>  \.J  -  S,  -  lU-=h,jlB  -  \,IU-VjId 


“  b2  f,  %  ”  Vj^D 

Hu  ^  )Ib  *-*ki  n  —  2  i 

l  '  1“  -  “h.jlB  J 


,  V  u  €  MU.  ). 

*1 


As  u  varies  over  MU.  )  it  is  clear  from  (3.30)  that 


»u  -  v,l 

h  -  Vjl 


minimized  for  that  un  that  minimizes  flu-u.  namely  for  u_  =  EU  )  u 

u  h ,  j "  d  u  k  i  n ,  j 

Thus  we  have 


-h 


(3.31a) 


ie(Wj-viId  .  h  -  vjip  5 

lE(Xk,>5h.j-=h.jlB  l“j-VjlB 


<  C^vi h). 


We,  of  course,  also  get 
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co  Q|c\i  cn 


(3.31b) 


Ie(V\j-\jIb 


<  C^th). 


Estimates  (3.31)  are  similar  to  (3.28),  but  involve  a  different  pairing  of 
approximate  and  exact  eigenvectors. 


Remarks  3.4.  Pierce  and  Varga  [12]  proved  eigenvector  estimates  in  ||*||  and 
BabuSka  and  Osborn  [6]  established  them  in  ||*||  . 
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L  Precise  Asymptotic  Estimates  for  the  Eigenvalue  and  Eigenvector  Error 

In  this  section  we  use  the  notation  introduced  in  Subsection  3.c),  i.e., 

we  let  k^  be  the  lowest  index  of  the  i^^  distinct  eigenvalue  of  (.2.5)  and 

suppose  A,  has  multiplicity  q  . 

M  A 

a)  The  Eigenvalue  Error 


For  i=l,2, •••  and  j=k^ , • • • , k^+q^-1  fixed,  Chatelin  [7,8]  has 
shown  that 


(4.  1) 


(A.  .-A.  )/A. 

h.  J  kt  kt 


H<I'Ph,E'Uk,,uh,jllB  /  ®E<Ak1,uh, jIb 


1,  as  h  0  (cf.  1.3c)). 


We  now  prove  a  refinement  of  (4.1)  (cf.(1.3c)  and  (1.4)) 


Theorem  4.1.  For  i  =  1,2, •••  there  is  a  constant 


(4.2) 


Uh,J  y^k! 


I(,-VeiVvjIb/  Ie'VVjIb 


-  ll  < 


dj.  such  that 
2 

d^  (h),  j=k. ,  •••,ki+qi-l, 


where  rj(h)  is  defined  in  (3.7). 


Proof  Let  u  *  E(A,  )u.  We  have 

-  k,  h.j 

(4-3)  (%-ph.j)B(u'uh.j)  '  B(Tu’Vj)  -  B<u’ThVj) 

=  B(u’(T-Th)Vj) 

=  B(T(I-P.)u,u)  *  B(T( I-P  )u.  u.  ,-u) 
h  h  h,  j 

=  B(T(I-Ph)2u,u)  +  B(T(I-Ph)u,uh  j~u) 
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=  B((I-Ph)u,  (I-Ph)u)  +  DUl-P^u^  j-u). 


Using  the  fact  that  B(u,  u^  =  B(u,E(A^  )u^  =  j|u{|g  ,  (4.3)  can  be 

written  as 

S7771  Ml  -  XT-  *  d<(i-Vu’Vj-u)- 

n,j  K*  K\ 

2 

Dividing  by  ||(I-Ph)u||g  ,  multiplying  by  *h  ,  and  subtracting  1  from  both 
sides  we  find 


(4.4) 


‘Vj'V'V  t  A  dii1'pI'I"A,j~“1 

l(I-ph,“l|/  Ml  x*.  h,J  |(>-ph)“»l 


From  (4.1)  or  the  standard,  well  known  results  for  eigenvalue  approxima¬ 
tion  we  have 

■\2 


(4.5) 


V  J~\  < 

*k.  1 


sup  1(I-P  )u|| 

w  /  -V  \  HD 


ueM(A,  ) 

Mb’*' 


=  d. 


^  sup  ||(I-Ph)Tu|| 


ueM(A.  ) 

IWIlX 


h  HB 


<  djl,  (h),  j  =  kj.—.kj  ♦  qrl, 


from  (3. 13a)  we  have 


(4.6) 


8"-ph)ulD  5  ,,(h)  »(1-ph,u»B. 


and  from  (3.26)  and  (3.31b)  we  have 


(4.7) 


I* 


h.j 


=  !uh,j"EUk1 


)Uh, jUd 
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S  dj  »(l>>IVj-EUk,)Uh.jllB 

-  d,  n(h)i(I-Ph)u|B  . 
Combining  (4.4)  -  (4.7)  we  obtain 


\jV\ 

|u-ph)ul|/  M| 


A  .  |D(  ( I-P.)u,  u,  -u| 

S  d,„2(h)  .  - h 


the  desired  result. 


<,  d^  + 


2  Ah, j  uHd 


H(I-ph)ullB 


*  djT ?, 


Remark  4.1.  Formula  (4.4)  Is  due  to  Chatelin  [5,6],  and  is  used  by  her  to 
prove  (4.1).  Using  eigenvector  estimates  in  ||’Hg  (0.26))  one  can  prove 


(Ah.j-y/Aki 


l(I-Ph,E(V'VjllB/  *E‘Xk.’VjlB 


-  1 


S  d^Qi). 


(4.2),  which  was  proved  using  eigenvector  estimates  in  |(*Id  (O- 31b)  together 

with  (3.26)),  is  an  improvement  over  this  result  since,  as  we  saw  in 
2 

Subsection  3.b),  n  may  be  of  higher  order  than  v. 

Theorem  4.1  relates  the  eigenvalue  error  (A,  -A.  )/A,  to  ||(I-P,  )u|| 

with  u  =  E(A^  )u^  j.  We  now  prove  a  result  that  relates  the  eigenvalue  error 

to  ||(I-Ph)u||g  /  (|u|||  ,  where  u  €  M(Aj^)  and  E^(A^ )  u  =  ,  i.e., 


u 


as  defined  in  Subsection  3.d). 
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C\J  03 


Theorem  4.2.  For  i  =  1,2, there  is  a  constant  such  that 


(4.8) 


tXh,J_Xki)/Xk1  ,  „  ...  .  « 

2  p  ~  1  ^  (h),  J  -  k. i i  k,  q,~l, 

I«-V“Ib/Wb  iii 


where  u  €  )  satisfies  )u  =  j- 


Proof.  With  u  e  M(X^  )  satisfying  E^(X^  )u  =  we  have 

“vvj’  B  (u'Vj’  ■  B(Tu’uh,j)  ■  B,uJh\j1 


=  B(T(I-Ph)u,u) 


+  B(T(I-Ph)u,uh  j-u) 


-  'Sclil!-vui! 


+  B(T(I-Ph)ulUh  j-u), 


from  which  we  get,  as  above, 


(4.9) 


Uh..1-V/Xlo 

l(I-ph)u«l  /  IEh(\1)u«i 


,  -  VjA,  ,  A  Dl(1-phln.j-ul 


h,J  |ti-ph)»i 


It  follows  from  (3.13a)  and  (3.16b)  that 


(4. 10) 


#Vj  ■  u»d  5  <yi(i-ph)uiiB 


Combining  (4.5),  (4.6),  (4.9),  and  (4.10)  we  obtain 


uh.,rV/xh, 

I(I-Ph)uiB  X  «Eh<:'k1,ulll 


<  di  rj  (h). 


from  which  we  get 
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<m  cq 


(4. 11) 


(x.  rK  )/x. 

h,  J  kt  kt 


|(I-Ph)u||  /  |u|| 


- — 2  5  diT)2(h) 


l»K 


■W^b 


Since  u  =  (u-E^(X^  )u)  +  E^(X^  )u  is  an  orthogonal  decomposition  in  Hg,  we 


have 


and  hence 


(4. 12) 


Hi*  i->yv»ii*  iEh(xit,,uiil 


Hb  .  ,  „  lu'E„<\,)“«i 


Using  (3. 16a)  and  (2. 13)  we  see  that 


(4. 13) 


I“-Eh(Xk, ,UIb  2 

-  -  i  “  <  C  cd(h) 

\t<L  U 


IEh(Xk,  )U«B 


S  C  i)  (h) . 


Combining  (4.11),  (4.12),  and  (4.13)  we  get  the  desired  result. 


b)  The  Eigenvector  Error 

Let  i  =  1,2, •••  and  let  j  =  -  1  be  fixed  and  consider 

j  and  E(X^)u^  ^  (recall  that  flu^  j||g  =  1).  We  showed  in  Subsection 
3.d)  (see  (3.31a))  that 


(4. 14) 


»E(xk,)5h,,fv  A £  di,,(h)#E(xk,)Vj-VjiiB 


From  Lemma  3.1  we  have 


(4. 15) 


Yj  -  xk,  ’  1eiV\j-VjIb  -  VPVvAjId 


Combining  (4. 14)  and  (4. 15)  we  obtain 


Xh.J  -  xk,  x  lE,Xk,)“h. j'YjIb  ,1-di’2(h»' 
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which  implies 


(4. 16) 


I^VVj-vA  . 


Xh,j  Xk, 


i(I-ph)EUk,)VjlB  ’  I(I-ph)E(Xkl)VjlB  (1-dl’>2) 


Since  j  =  E(\k  j  +  j“E(X^  )u^  j)  is  an  orthogonal  decomposition 
in  Hp,  we  have 

1  =  lE(Xk,>VjlD  *  K.J-  E(Xk,>Vj»D- 


From  this,  (3.26),  and  (4.14)  we  get 


(4.17)  1  s  |E(*  ky\jl  *  dl’>2|(Ph‘I)E(Xk1)"h,jlB  -  (1V> 

M 

Now,  combining  (4.2),  (4.16),  and  (4.17)  we  have 


ii^Vvj-VjIb 

*(I-Ph)E(VVjlB 


(xh..rxk,1/xk, 


1+d.Tj" 


l(I-ph,E(VYj»l  /I«VVj£  i-v2J 


1/2 


<  i 


,  dV>',/2 

(1  *  V  >  TT-2~ 


1-d^ 


.  ri  ,  2,1/2 

5  [1  +  CL  7)  ] 


<  1  +  djij  (h). 


We  summarize  this  (cf.  (1.3b)  and  (1.4))  in 


Theorem  4.3.  For  i  =  1,2, •••  there  is  a  constant  d^  such  that 


(4. 18) 


||E( Xr  ) j-uh  jllg  2 

!  *  - -1— 2  1  ♦  d/(h),  j  =  k  ,...,k  +q  -1. 

l(I-ph)E'Xkl)dh.jlB 


j  can  be  replaced  by  ^  in  (4.18). 
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Remark  4.2  (4.18)  is  stronger  than  (3.26)  since  y  may  be  of  higher  order 
than  v. 


Next  consider 
,-h 


u.  .  and  (recall  that  u*)  €  M(X,  )  satisfies 
n,  J  j  J  kj 


)Uj  =  j).  We  know  (see  (3.28a))  that 


,-h 


-h 


•“j  -  Vj«D  *  V!h>  »UJ  -  Vj»b 


This,  together  with  Lemma  3.1,  yields 


which  implies 


\,j  -  v  -  i=5  -  vj»b  -  \,i=5  - 


2 1=5  -  v/b  “-v2>- 


(4. 19) 


Huj  “h,A  s _ 

ja-ph)uj8|  (1'di7,2) 


Uh,J~  Akt)/Akt 

»(I-V-J«2B  „2  , 

- v— 5 -  ( 1-d, y  ) 

»“5«b 


Finally,  combining  (4.8)  and  (4.19)  we  have 

,-h  - 


,  s  lYVja 


1+d^' 


l-djU' 


1/2 


5  l  +  din. 


This  result  (cf.  (1.3a)  and  (1.4))  and  the  related  result  (3.19)  are 


summarized  in 


Theorem  4.4. 


(4.20) 


For  i  =  1,2, *••  there  is  a  constent  such  that 

l“r“h  Jb  iVEh(xk,  2 

1  £  — P’l .u  =  — J — sL"  <  1  +  d,7)  (h), 

l(I-Ph)53«B  i(I-Ph>=jiB 


J  =  W'- 

(4.20)  remains  valid  if  ^  is  replaced  by  u^  j  and  u||  by  u^. 
a  constant  d^  such  that 


There  is 


(4.21) 


Hu_Ehuk1  )uHb 

i  <  - ^ - -  £  1  +  d, u(h) , 

•(“-Ph>“llB 


for  all  u  e  M(A.  ). 

Ki 


Remark  4.3.  We  have  restated  (3.19)  in  (4.21)  because  it  is  related  to  (4.20) 
and  it  is  the  strongest  known  result  of  its  specific  type.  It  should  be  noted 
that  (4.21)  is  true  for  all  u  e  M(X  ),  whereas  (4.20)  is  valid  only  for 
u  =  u||,  j  =  kj.^’.kj  +  q^-1.  However,  for  these  u’s,  (4.20)  is  stronger 
than  (4.21). 


Remark  4. 4  See  [2,4]  for  a  numerical  study  of  the  reliability  of  the  results 
of  this  section  -  which  are  of  am  asymptotic  nature  -  as  a  guide  to  practical 
computations  -  which  often  takes  place  in  the  preasymptotic  phase. 

5.  An  Additional  Result  for  Multiple  Eigenvalues 
Theorem  4. 1  shows  that 


<  C 


inf  |EU  )Uh 
*€Sh 


ki 


iE'VYk' 


2 

B 


In  [3]  BabuSka  and  Osborn  proved  the  stronger  result  (cf.  (1.5)) 
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h.ki 


Xk,  5  C 


inf  inf  |u-Z|; 


ueM<V  Z€Sh 


IuIb=1 


(as  well  as  similar  estimates  for  A.  .  -  A.  ,  J  =  k.  +  l,***,k.  +  q  -1,  and 

n,  j  Kj  l  li 

for  the  eigenvector  errors),  which  shows  that  A.  .  -  A.  ,  the  error  in  the 

h,  kj  Ki 

approximate  eigenvalue  closest  to  A,  ,  is  governed  by  the  approximabi 1 ity  of 
the  exact  eigenvector  corresponding  to  A.  that  can  be  best  approximated  by 

M 

S^.  In  this  section  we  give  a  simplified  proof  of  the  results  of  [3],  which 
in  addition  provides  information  on  C  (the  results  in  [3]  only  established 
that  C  is  a  constant),  and  estimate  the  eigenvector  error  in  ||*||q  and 

IHI-b- 

As  above,  for  i  =  1,2,  •••  suppose  k^  is  the  lowest  index  of  the  i*"*1 
distinct  eigenvalue  of  (2.5)  and  let  be  its  multiplicity,  i.e.,  suppose 


X-i^i-r1  Aki_1  <  Aki  =  Ak,+1  ","*Xk»+vl<  Xki+cn c  xki*i' 


Let 


(5.  1) 


e  .(h)  =  inf  inf  ||u-*|| 
U€M(A,_  )  ° 


k,'  *€Sh 


l|u«B=1 


B(u-Uh,k1)s,,,=B(u'\k1+j-2)  = 


=  0 


inf 

ueM(A.  ) 

Me"1 


eu(h) 


^“•Vki1 . B(u’\k1+j-2>  =  °* 


* '  qi  ’ 


where  M(A^)  is  defined  in  (2.9).  The  restrictions  B(u,u^  k  ) 

=•••=  B(u,u^  =  0  are  considered  vacous  if  j  =  1.  We  note  that  they 

are  equivalent  to 


B(u, E(Ak^)uh  ^)  =  0,  +  j  -  2 
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and  to 


B(Eh‘ V“'Vt>  *  *  j  - 2- 


Theorem  5. 1  (cf.  (1.5))  For  i  =  1,2, •••  there  is  a  function  C^(h)  and  a 


constant  C^,  with 


(5.2) 
such  that 


C^(h)  i  1  +  d^t/(h),  d^  =  constant. 


(5.3) 


(K 


, ki+J-1  ~  Xkt)/Xkt  *  ^i(h)  ei,J(h)’  J*1*  “‘•‘V 


and  such  that  the  eigenvectors  ui*u2’“*  (2.5)  can  be  chosen  so  that 

(2.6)  is  satisfied  and  such  that 


(5.4) 


HV^+J-1  '  5  ^i(h)  ei,J(h)’  J=1****'c?i1 


(5.5a) 

and 

(5.5b) 


HVki+J-1  “  “ki+j-l^D  *  CiT,(h)ci,J(h)’  J  1,‘**,qi’ 


H^.ki+j-l  “  ^j+J-J-B  5  Ci*,(h)ci, j(h)» 


where  u(h)  and  n(h)  are  defined  in  (3.7)  and  (3.8). 

Proof.  Let  i  and  j,  with  i  =  1,2, •••  and  j  =  l,***,q^  be  fixed.  Note 
that  e  (h)  £  A  v(h),  Vue  M(A.  )  and  e,  .(h)  S  A.  vfh),  j  =  l,***,q.. 

u  Kj  Kj  1 1  J  Kj  1 

Let  ueM(Aki)  with  B(u,uh  ki)  =  •••  =  B^.u^^)  =  0  and  |u|B  =  1. 
Now  apply  (2.7*1)  emd  Lemma  3.1  with  (A 

Eh(Xk,)u 


■u)  *  (v 


sind 


w  = 


Since 


by  the  orthogonality  of  the  approximate  eigenvectors,  and 
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B(E^(X^  )u, t)  -  B(u, |)  =  0,  l  =  ***»ki+j-2 


by  the  assumption  on  u,  we  have 


(s-6)  Wj-i  '  xki  s  b 


r  W"  w* 


IiW-Id  IVV-IdJ 


-  *. 


|EhUk,)u-u||-Ak,»Eh(Ak,)u-ulg 

lEh(Xk1,ullD 

,  l^^k,  )u-»l 

«VXk,)U«D 

From  (3. 19)  we  have 

(S.7)  I^'V^Ib  S  ( 1+dw) |u-Phu|g  . 

From  (3. 13)  and  (3. 16b)  we  see  that 


I  IVV-Id  -  xr>  '  HW’-lp  -  MdI 

5  lEh(XkI,u-“lD 

5  0  «u  -  V»D 

s  dn  ||u-Phu|lB 
=  d7j(h)  £u(h). 

which  shows  that 


(5.8) 


Nc^'V^D 


1  + 


drje^h) 


-»2 


lEh(\1)ullDJ 


i  1  +  di?  eu(h) 
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£  1  +  di)  v. 


Combining  (5.6)  -  (5.8)  we  get 

(5'9)  5  11  *  dr)zU*dm)  |u-Phu(|| 

s  <1  *  <*>  |u-VIb  ■ 


Now  since  (5.9)  holds  for  all  u  €  M(A^  )  with  EKu.u^  ^)  =  0, 


l  =  ••*,ki+j-2,  and  |ju||B  =  1,  we  have 


'Vk^j-rV^k, 4  (Udv(M) 


inf  Inf  ||u-*|| 

ueM‘V  «Sh 

Mb-1 


B 


B(u,uh.k,) . !,u'\vj-2W 


=  ( 1+d^yfh)  )e^  j(h). 


which  is  (5.3)  with  C^(h)  =  l+d^u(h).  Thus  (5.2)  and  (5.3)  have  both  been 
proved. 


Remark  5.1.  (2.7**)  and  Lemma  3.1  lead  to  a  particularly  simple  proof  of  a 

result  slightly  weaker  them  (5.3)  for  the  case  i  ■  J  ■  1.  It  follows 
immediately  from  these  two  results  that 


Xh,l  ~  X1  *  B 


P.u  ’ 
n 


V 


-  X. 


«phu-uii|  - 

»Phu»D 


iPhu'u#B 

S  — - -  V  u  €  MUt), 

«Ph“lD  1 
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and  hence 


(A.  ,-A.l/A.  $  inf 
n,  l  l  l 


URu-ul 


U€M(A.)  A.||P.u||~ 

8u«B=1 


|  <  C(h)  e1  jth). 


where  C(h)  — »  1. 


Now  consider  (5.4)  and  (5.5).  Let  i  =  1,2, •••  and  j=l , - • • , be 
fixed.  Let  u^j.,  e  MU^  >  satisfy  VVWl  =  “h, k.tj-l' “i,t j-1  = 
where  +^_1  was  introduced  in  Subsection  3.d)).  Applying  Lemma 

3. 1  with 


we  get 


(A, u)  = 


Vj-1 


[ ki’  Kk^-ih 


amd  w  = 


Vki+J-1 
8^,  kt  +j-i  8d 


(5. 10) 


h> \  -  II  ra  i\kl;HiD  «b 


ki  Kk^j-ifc 


Vh  If  . 

^.ki+J-l^D  ||D 


From  (3.28b)  we  have 

(5.H)  I'Scj+j-!  "  uh, ki+j-1  Hd  S  CT,(h)K1+j-l  ‘  Uh, ki+ j~l Hb  ' 


(5.  10)  and  (5. 11)  yield 


A 

h, ki+j-1 


A 

ki+J-l 


> 


[1-A.  C27j2(h)] 

_ M _ 

llVki+j-Jo 


K.ki+j-i 


which,  together  with  (2.11),  (2.12),  and  (5.3),  yields 


37 


(5. 12) 


K.+J-l  "  Uh, ki+J-1  Hb  5 


(X  -  X  )1/2 

Uh,kt+J-1  AktJ 

.  1/2  r.  .  _2  2,1/2 

Xh.k,*j-1  "V" 


X*/2C  (h)1/2e  .(h) 

*i  1 _ 1 .  J _ 

.1/2  ri  _2  2,1/2 

V»,.j-i  u'\,c  ” 1 


<  C^h)  e  (h),  j  =  1, 


.q 


i’ 


where,  because  of  (5.2),  C,(h)  <  1  +  d i>(h).  (5.12)  shows  that  the  u'  .  , 

i  i  kj+j-l 

satisfy  estimates  (5.4).  (5.12),  together  with  (3.28b)  and  (3.29b),  shows 

that  the  +j_^  satisfy  estimates  (5.5).  They  will  not  in  general, 
however,  be  orthonormal  with  respect  to  B,  so  that  (2.6)  may  not  be 
satisfied. 

It  remains  to  modify  the  u/  ,  i.e.,  replace  u'  by  u,  , 

Ki+j-i  k1+J-l  m+j-1 


in  such  a  way  that  (2.6)  and  (5.4)  and  (5.5)  hold.  We 

“i, 


j.  Let  j  =  1.  If  we  define  u,  =  n-,— n  , 

11^ 

(2.6)  is  satisfied.  From  (5.12)  we  have 


we  have 


proceed  by  induction  on 

|\  |B  “  !»  so  that 


(5131  *  |[>  *  -  Vk,.j-1«B]1/2-  ‘I 

<  ~  Uh, kj+J-1  I'b 

2 

<  Cu  ei  J(  J-l.—.qlf 

and  hence 


Uh, kt  Hb  ~ 


u. 


h,  k, 11 B 
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s  IK,  1b'1  I  *  K, '  \h,Ib 

5  Cyci , i 1  £i(h),,ci,i 

£  Ci(h)ei  1Ch), 

where  C^(h)  £  1  +  d^v(h),  which  is  (5.4)  for  j  =  1.  Using  (3.9),  (5.13), 
and  the  fact  that  the  u,'  ,  .  satisfy  (5.5a)  we  get 

Ki+j-1 

K,  '  Vic.IId  £  K, '  “i.lD  *  K, '  “h.k.lD 

■^1  K,«b'1I  *  K,'\k,l» 

s  Cpci,i  *  Cwi,i 

i  CjTjthiEj  j(h), 


which  is  (5.5a)  for  j  =  1.  A  similar  estimate  establishes  (5.56)  for  J  *  1. 
Next  suppose  j  ~  2.  Let  u£  =  uj^  +  1  “  B(u^  +  ^ ,  )u^ .  Using  (5.4) 


for  J  =  1,  (5.5b)  for 

u,'  ,  ,  and  that  e.  . 

Ki+j-1  i, 


j  =  1,  (5.12),  and  the  facts  that  (5.5b)  holds  for  the 
^  e.  ,  we  have 

1  y  ^ 


(5.14)  |B(uici,1.uki)|  <  *  ISIVlc.rVrWk.’l 


*  iBtuk,*i'ukrVk,i 


=  xk,lDK,.rVk,.i’uk1,l 


*  lBluh,k,*ruk|*l'ukruh,ki 1 1 


*  Vi|mVi,ukfuh,k1)l 
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5  KcJb 

+  iVkt+r^+JB  IWkjB 
+  Xk1+l®Ukl  +  lHB  H'WkJ-B 

S  XkJuk1+ruh,ki+J-B 

+  K,k1+i  -  Vik  KrVkjB 
+  Xkt  (1+Buk1+l-uh,kl*JlB)lluk1-uh,k1ll-B 

SC^i,2  +  Ci(h)ci2Ci(h)ci(1  ♦C*C1(h)e1>1 
<  Cw(h)ci  2(h), 

and  hence 

(5.15)  Kitl-u{iM|B  =  S  Ct»(h)Cj  gth) • 

ll£ 

Now  set  u^+j  =  jj^  1 — |j— .  Combining  (5.12),  (5.13),  and  (5.15)  we  obtain 
®\1*ruh,ki*l®B  5  K^iKimIb  *  K.i'VmiIb 

*  I  *  K,.r  Vic,.Jb 

s  IK, .,1b-1  I  *  2K,.,-  \JB 

*  K,.r  Vk,.,lB 

*  C^Cj  2(h)  ■»  Cl,  Cj  2(h)  ♦  CjlhlCj  2(h) 
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*  C^hle,  2(h). 

where  C^h)  £  1  +  d^iKh),  which  is  (5.4)  for  J  =  2. 

Now  consider  (5.5a)  for  J  =  2.  Using  (5.13),  (5.14),  (5.15)  and  the 
fact  that  the  u£  satisfy  (5.5)  we  have 


l^i+rVki+JD  "  luk|+l_uk1+llD  +  l^+l  Vki+lt 

=  \\/Z  I  1%+Jb"1!  +  Kj  +  rVk^llD 


s  \\/2  I  luki+lb“1l  +  \\  l^+l’^+lb 
+  K,+1  -  VmiId  + 


+  Cyel,2  +  Cpe i , 2  +  Cvcl,2  +  Cyci,2 


5  Clvcl2( h), 

which  is  (5.5a)  for  J  =  2.  The  proof  of  (5.5b)  is  similar. 

Continuing  in  this  manner  we  get  (2.6),  (5.4),  and  (5.5)  for  j=l,*»*,q.. 
This  completes  the  proof.  o 


Theorem  5.2.  For  i  =  1,2,  •••  there  is  a  function  C^h)  with 

(5.16)  C^h)  >  1  -  d^th),  di  >  0,  constant, 

such  that 


(5. 17) 
and 

(5. 18) 


Uh,k!+j-l  '  \]/\  '  Ci(h)  Ci , J(h) ’  J  =  1’***’qi* 


H^.kt+j-l  "  Ukj  J— 1  "  S(h)  Gi ,  j(h^  ’  J  ~  1,*‘*,qi 


Proof.  First  consider  (5.18)  for  J  *  1.  It  is  immediate  that 
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2:  inf  inf  (ju~zj 


U€M(X.  ) 
*1 

luiB=1 


z«s. 


B 


■  ct,i'h)- 


Thus  for  J  =  1,  (5.18)  holds  with  C^h)  =  1. 

Now  suppose  J  =  2.  Since 

B(uki+l,uh,k1 }  =  B(uk,+l,Eh(;Sci  )uh,k,) 
•B(Eh\)ViiUh,k1) 


=  B(uh.k1+l*uh,k,) 
=  0, 


we  see  that 


*WrVtlB  1  Utiu.  > 

Mb*1' 

B(u,uh.k,)-0 


inf 

xeSK 


MR  = 


ei>2(h) 


Combining  this  result  with  (5. 13)  and  (5. 15)  we  get 

^uh,kl+ruk!+ Jb  “  Huh,k1+rukl+ Jb  "  ®uk|+l“uk1+l Bb  "  luk1+l“uk1+llB 

“  luh,k|+l“uk1+llB  “  2^uk1+ruk1+l^B  “ 

£  (l-diy)ci  2^h) , 

which  is  (5.18)  for  j  =  2.  Continuing  in  this  manner  we  get  (5.18)  for 
J-l.~-.qj. 

Now  consider  (5.17).  From  Lemma  3.1,  (5.5a),  and  (5.18)  we  see  that 


Uh,kl+j-l 


=  l^kt+J-l  ~  %+J-Jb 

*kJVki+J- l^D 


IK.kt+J-l  -  “ki+J-Jp 
BVkj+j-Jo 


>  h,kt+J- 

"  V 


1  [(i-di*'>2  -  ■ 


which  implies  (5.17). 


Remark  5.2.  Note  that  in  Theorems  5.1  and  5.2  we  have  shown  that 


Uh,k1+j-l  Xk,)/Ak, 


C'j(h) 


-  1 


£  dj,i>(h), 


whereas  in  Theorems  4.1  and  4.2  we  showed  that 


(Yj  -  V7Y 


«(,-Ph,E(:'1c1)Yj»B/IE(Y)Yj«B 


-  1 


<  diY)  (h), 


and 


(Yj  -  V7Y 


-  1 


S  dji)  (h). 


l<I-Ph)u|E/|u|| 


for  u«Muki)  with  yyu.Vj. 


Remarks  5.3.  For  a  computational  illustration  of  the  results  in  this  section 
see  [3,4]. 
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The  Laboratory  for  Numerical  analysis  is  an  integral  part  of  the 
Institute  for  Physical  Science  and  Technology  of  the  University  of  Maryland, 
under  the  general  administration  of  the  Director,  Institute  for  Physical 
Science  and  Technology.  It  has  the  following  goals: 

o  To  conduct  research  in  the  mathematical  theory  and  computational 

implementation  of  numerical  analysis  and  related  topics,  with  emphasis 
on  the  numerical  treatment  of  linear  and  nonlinear  differential  equa¬ 
tions  and  problems  in  linear  and  nonlinear  algebra. 

o  To  help  bridge  gaps  between  computational  directions  in  engineering, 
physics,  etc.,  and  those  in  the  mathematical  community. 

o  To  provide  a  limited  consulting  service  in  all  areas  of  numerical 
mathematics  to  the  University  as  a  whole,  and  also  to  government 
agencies  and  industries  in  the  State  of  Maryland  and  the  Washington 
Metropolitan  area. 

o  To  assist  with  the  education  of  numerical  analysts,  especially  at  the 
postdoctoral  level,  in  conjunction  with  the  Interdiscipl inary  Applied 
Mathematics  Program  and  the  programs  of  the  Mathematics  and  Computer 
Science  Departments.  This  includes  active  collaboration  with  govern¬ 
ment  agencies  such  as  the  National  Bureau  of  Standards. 

o  To  be  an  international  center  of  study  and  research  for  foreign 

students  in  numerical  mathematics  who  are  supported  by  foreign  govern¬ 
ments  or  exchange  agencies  (Fulbright,  etc.) 

Further  information  may  be  obtained  from  Professor  I.  Babuska,  Chairman, 
Laboratory  for  Numerical  Analysis,  Institute  for  Physical  Science  and 
Technology,  University  of  Maryland,  College  Park,  Maryland  207^2. 


